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Abstract 

A continuous measurement of energy which is sharp (perfect) leads 
to the quantum Zeno effect (freezing of the state). Only if the quan- 
tum measurement is fuzzy, continuous monitoring gives a readout E(t) 
from which information about the dynamical development of the state 
vector of the system may be obtained in certain cases. This is studied 
in detail. Fuzziness is thereby introduced with the help of restricted 
path integrals equivalent to non-Hermitian Hamiltonians. For an oth- 
erwise undisturbed multilevel system it is shown that this measure- 
ment represents a model of decoherence. If it lasts long enough, the 
measurement readout discriminates between the energy levels and the 
von Neumann state reduction is obtained. For a two-level system un- 
der resonance influence (which undergoes in absence of measurement 
Rabi oscillations between the levels) different regimes of measurement 
are specified depending on its duration and fuzziness: 1) the Zeno 
regime where the measurement results in a freezing of the transitions 
between the levels and 2) the Rabi regime when the transitions main- 
tain. It is shown that in the Rabi regime at the border to the Zeno 
regime a correlation exists between the time dependent measurement 
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readout and the modified Rabi oscillations of the state of the mea- 
sured system. Possible realizations of continuous fuzzy measurements 
of energy are sketched. 

PACS: 03.65 Bz 



1 Introduction 

Rabi flopping is a well known phenomenon: If a two-level atom with energy 
levels E\ and E% is under the action of a monochromatic electromagnetic field, 
which we assume for simplicity to be resonant, then the state vector of the 
atom oscillates between the lower and upper state with the Rabi frequency 
which depends on the amplitude of the driving field. Speaking more precisely, 
the probability for the atom to be, say, on the level E\ oscillates sinusoidally 
between 1 and 0. Whenever a von Neumann ("instantaneous" and perfect) 
measurement of energy is made at a certain time, the measurement result 
is either Ei or E 2 with the probability depending on the time when the 
measurement is performed, in accordance with this sinusoid. The atom is 
after the measurement in the energy eigenstate corresponding to E\ or E 2 - 
This means that such a measurement disturbs the Rabi oscillation radically: 
the phase of the Rabi flopping during next period is determined by the result 
of the von Neumann measurement. 

It is also well known that if perfect measurements of energy of the atom 
(von Neumann measurements) are very frequent and in the limit continuous, 
the interaction with the measuring apparatus prevents atomic transitions 
between levels. The atom stays in an eigenstate and the readout of the 
continuous perfect measurement reduces to the constant value E x or E 2 . 
This is the often discussed quantum Zeno effect |I| which has been verified 
experimentally 0. It is impossible in this regime of measurement to obtain 
information about the Rabi oscillation. Moreover, the Rabi oscillation is 
prevented by the measurement because the influence of the measurement on 
the atom is stronger than the influence of the driving field. 

In contrast to this we want to address in the following the fundamental 
question if it is nevertheless possible to make the Rabi oscillation "visible" 
by continuous measurement of the atomic energy, i.e. to obtain a correlation 
between the measurement readout E(t) and the Rabi oscillation of the state 
vector. To achieve our goal we must first of all avoid the Zeno effect by 



2 



weakening the influence of the measurement. This can be done with the 
help of continuous measurements which are fuzzy (unsharp) instead of being 
perfect (sharp). 

A fuzzy (unsharp) measurement is caused by the low resolution of the 
respective detector. In the case of an observable with a discrete spectrum 
fuzziness means that the measurement is not able to transfer the state of the 
system to an eigenvector. Fuzziness does not go back to a defect detector, 
but to its finite resolution. Note that in practice a fuzzy measurement is 
the generic case and that a sharp or perfect measurement may be considered 
only as an idealization. The study of continuous fuzzy measurements opens a 
very interesting new field in which fundamental problems of quantum theory 
as well as potentially important concrete applications are discussed. There 
are physical situations in which we have to rely on fuzzy measurements if we 
want to obtain relevant information M. 

The theory of continuous quantum measurements has a long history. An 
incomplete list of references is ||, f| ^| 0, ||, [|, 0, [0], [12|, [13], [14j| . The majority 
of authors studying continuous quantum measurements base their investiga- 
tions on particular micro-models for the measuring device, see || || |K| |J for 
the measurement of the energy of a two-level system. As compared to this 
the phenomenological and therefore model-independent approach presented 
below has the big advantage that the universal structures of the theory of 
continuous measurements can be revealed, including those connected with 
fuzziness. This approach based on restricted path integrals (RPI). It has 
been proposed in || [HJ. The approach has already been applied to the mea- 
surement of energy in |16 , \L% . However, an ad hoc assumption was made in 



n| [17| that the measurement readout is a constant curve E{t) = E n = const. 



This is too restrictive and misleading for a continuous fuzzy measurement. 

A theory of continuous fuzzy measurements (say, monitoring some observ- 
able has to answer the following questions: How can the type and strength 
of fuzziness be represented theoretically and what will be the character- 
istic parameters? Let the measurement of the observable A(p,q,t) be of 
duration T, what is the probability to obtain the measurement readout 
[a] = {a(t)\0 < t < T}? If for an initial state \i/jq) of a system a partic- 
ular measurement readout [a] is obtained, what will be the state \ip^) °f the 
system at the end of the measurement? The way to answer these questions 
in the framework of the phenomenological RPI approach will be presented 
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in Sect. |2|. The resulting dynamics is governed by a complex Hamiltonian 
depending on [a]. The probability for [a] can be read off from |V>t')- We 
obtain as a characteristic property of a continuous fuzzy measurement that 
its resolution improves with increasing duration T as Aa^ ~ 1/ yT. This 
reflects a Zeno-type influence which is still present even if it does not result 
in the proper Zeno effect 

In Sect. |H the scheme is applied to the measurement of energy of a mul- 
tilevel system which is not influenced otherwise (no additional driving field). 
Different durations T and the corresponding regimes of measurement are 
discussed. The time after which neighboring levels are resolved is called the 
level resolution time T\ r . Its shortness is a measure for the strength of a 
Zeno-type influence ("level freezing") present in the continuous fuzzy mea- 
surement. The relation to the process of decoherence leading to the perfect 
(von Neumann) energy measurement is studied. 

A two-level system under the resonant influence of a periodical driving 
field is then investigated in detail in Sect. f|. The driving field acts in oppo- 
sition to the Zeno-type influence of the measurement lasting for the time T. 
Short Rabi periods Tr ~(Rabi frequency) -1 represent strong influences of 
the driving field. It is obvious that different relative magnitudes of the times 
T\ T , Tr and T characterize different regimes of measurement with different 
behavior of the atom state. A complex general view is to be expected. We 
show that there is a regime in which a correlation exists between the time 
dependent measurement readout E(t) and the Rabi oscillations of the atom 
state under the influence of the periodic driving field. 

In Sect, [5] possible schemes of realization of the measurement in question 
are discussed. A continuous fuzzy measurement may be realized, according 
to these schemes, as a series of "instantaneous" (short-time) unsharp mea- 
surements. Finally in Sect. |6] concluding remarks are made. 

2 A phenomenological approach to continu- 
ous fuzzy measurements 

1 This phenomenon is referred to as "narrowing" by Milburn in the second of Ref. H. 
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2.1 Restricted path integrals 

Evolution of a quantum system when no measurement is performed (and 
therefore the system is closed) can be described by a unitary evolution op- 
erator U t . The matrix element (q"\U t \q') = U t (q",q') of the evolution op- 
erator between states with definite coordinates (called the propagator) may 
be constructed with the help of the Feynman path integral (for example in 
the phase-space representation, i.e. with integration over paths \p, q] in the 
phase space). 

The formalism of path integrals is physically transparent and can be 
generalized to different non-standard situations. One of them is the situa- 
tion when the system undergoes a continuous measurement (and is therefore 
open).[| The idea to apply restricted path integrals (RPI) in this case was 
suggested by R. Feynman JTJ) and technically elaborated by one of the present 



authors (see |J for the original papers and |15j for modern formulations). 



A continuous measurement performed in the time interval [0, T] gives 
some information about which paths the system propagates along. The cen- 
tral idea of the RPI approach is that the original path integral must be 
changed in accordance with this information. In the simplest case the path 
integral has to be restricted to a subset of paths depending on the mea- 
surement readout. In a more general case the contributing paths should be 
weighted by a weight functional depending on this readout. 

To be more concrete, let us consider as a typical continuous fuzzy mea- 
surement the monitoring of an observable A(p, q, t) during the time interval 
[0,T]. In this case the registered measurement readout (output) is a curve 

[a] = {a(t)\0 <t< T}. (1) 
Then the path integral has the form 

U?(q",q') = J d[p]d[q}w [a] [p,q}erJo^~ H ^\ (2) 

The weight functional W[ a ] will be discussed below. 

Evolution of the system undergoing continuous fuzzy measurement re- 
sulting in a certain readout [a] may now be represented by the following 

2 Of course, the measurement must be precise enough so that quantum effects be essen- 
tial; then it is called a quantum measurement. 
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formulas: 

l#> = t4%o>, pP = U?p (u^ ] )\ (3) 

These formulas contain partial propagators depending on [a]. The first of 
these formulas may be applied to the evolution of pure states while the second 
is valid both for pure and mixed states. Note however that both of them are 
applicable only in the case when the measurement readout is known, i.e. in 
the situation of a selective measurement. A pure state remains pure after the 
selective measurement. 

The norm of the state resulting from the measurement is, according to 
these formulas, less than unity. Normalization may of course be performed, 
but it is much more convenient and transparent to work - as we will do - with 
non-normalized states. Their norms give according to (jTSJ the probability 



densities P[a) = TrpJp of the corresponding measurement readouts [a]. The 
following formula (with an appropriate measure of integration): 

Prob([a] eA)= J d[a] P[a] = d[a] Tr (u?p ■ ( 4 ) 

gives the probability for a measurement to result in an readout from the set 
A. The word 'density' refers to the space of all measurement readouts. 

In the considerations above we discussed the situation when the mea- 
surement readout is known (selective description of the measurement). If it 
is only known that the continuous fuzzy measurement is performed but its 
concrete readout is unknown (for example in a priori calculations), we have 
the situation of a non-selective measurement. In this case the evolution is 
described by the second of the formulas @ with summation over all possible 
measurement readouts: 



Pt 



d[a)pP= d[a)U? Po (u?y. (5) 



The density matrix px is normalized for an arbitrary initial (normalized) 
state po if the generalized unitarity condition 

J d[a] (U [ * ] f 4 a] = 1 (6) 

is valid. 
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In the most general case we have only a partial information about the 
measurement readout. If we know that [a] G A, the evolution of the system 
must be described by the formula 

p$ = J a d[a) pP = J a d[a] U? Po (U?y (7) 

intermediate between Eq. (Q) and Eq. (|5|). 

To specify the type of fuzziness we are dealing with, we have to specify 
the weight functional W[ a ] in (0). To do so, we introduce the concept of 
closeness in the space of curves [a] characterizing readouts of the continuous 
fuzzy measurement. By being close we shall mean that the temporal mean 
squared deflection abbreviated according to 

((A - af) T = I £ (A(p(t), q(t),t) - a(t)) 2 dt (8) 

is small. According to the Feynman approach A (p, q, t) is a c-number func- 
tion corresponding to the quantum observable. @ is an integral characteris- 
tic of the difference between two functions. This characteristic is large (and 
the curves are therefore not close) only if the deflection of one curve from 
another is big during a long time. The curves which have large deflections 
but only during short intervals are considered to be close. Q 

Based on this, the information supplied by the readout [a] of a contin- 
uous fuzzy measurement can be typically presented by the Gaussian weight 
functional 

W[a) [p, q] = exp 

which is exponentially small if ((A— o) 2 )t is big. The inverse of the parameter 
K may serve as a measure of fuzziness. Later on, it will be connected with 
the resolution of these devices. Eq. (||) characterizes a particular class of 
measuring devices.^ 

3 We accept here all curves [a] as possible measurement readouts. In fact only com- 
paratively slowly varying curves [a] must be considered if the inertial properties of a real 
measuring device are taken into account. Then large deflections for a short time are 
impossible. 

4 The case of a sharp continuous measurement (when a single path presents a measure- 
ment readout) is considered in the paper of Aharonov and Vardi ]l9| . 



(A(p(t),q(t),t)-a{tj)' 



dt 



(9) 
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We shall consider the case when the fuzziness of the measurement does 
not alter with time, so that k is constant. In principle, the parameter k 
could depend on time (then it must be included in the integrand of the time 
integral). In any case k cannot depend on the duration T of the measurement. 
This follows from group-theoretical properties of continuous measurements 
[|T5| , |20| , but can also be derived from models of measurements [14]]. An 



important consequence of this fact will be discussed below, namely that the 
resolution of the continuous measurement is improved according to Aay ~ 
1 / yT when its duration increases. 



2.2 Complex Hamiltonian 

Using the weight functional (H) in Eq. (§), we find for the partial propagators 
describing the evolution corresponding to a certain measurement output [a] 

U? (<f, q') = J d[p] d[q] exp j % - £ (pq - H [a] (p, q, t)) dt\ (10) 

where p = p(t) and q = q(t). The restricted path integral (|T0| ) has the form 
of a non-restricted Feynman path integral but with an effective Hamiltonian 

H [a] (p, q, t) = H{p, q, t) - ikTi {A{p, q, t) - a(t)f (11) 

containing an imaginary term.Q 

The possibility to describe the dynamics by an effective Hamiltonian is 
very important since it makes practical applications rather simple: one may 
forget about path integrals and reduce the problem to the solution of the 
equivalent Schrodinger equation with a complex Hamiltonian. It is impor- 
tant also from the theoretical point of view because the imaginary part of 
the Hamiltonian indicates a characteristic trait of the process of continuous 
measurement namely that information is dissipated.^ 

5 The form of the imaginary term may be different for another weight functional wr a i , 
i.e. for another class of measuring devices. 

6 Let us remark that restricted path integrals describing continuous fuzzy measurements 
and leading to non-Hcrmitian Hamiltonians must not be confused with path integration 
in the presence of infinite potential barriers. Potential barriers also do not permit paths 
to go out of some corridor, but the unitary character of the process is not violated, no 
dissipation occurs. 
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2.3 Measurement of energy 

We turn to the measurement of energy. Let the Hamiltonian H of the system 
have the form H = H + V where H is the Hamiltonian of the "free" 
multilevel system and V is a potential describing an external influence leading 
to transitions between levels if no measurement is applied. By measurement 
of the energy we shall mean the measurement of the observable A = Hq 
leading to the readout [E] = {E(t')\0 < t' < T}. 

The weight factor (|^) with A = Hq and [a] = [E] can be written in the 
form 



w \E}\p,q\ = exp 



((H -E) 2 ) T 
AE 2 T 



(12) 



where we made use of the notation (|J) and introduced the parameter 

AE T = — L= (13) 



characterizing the fuzziness of the measurement in a more transparent way 
than k. 

The weight factor (|T2"D guarantees that only those Feynman paths [p, q] 
are taken into account, for which the time dependence of the observable Ho 
is presented by a curve close to the curve [E] . The measure of the deflection 
((Ho — E) 2 )t is, according to (|12D , less or of the order of AE 2 ^.. This means 
that the resolution of the measuring device described by our formulas is equal 
to AE T .f\ 

It has been stated above that the parameter k does not depend on the 
duration T of the measurement. Therefore, it follows from Eq. (|T3| ) that the 
resolution of the continuous fuzzy measurement is improved when the dura- 
tion of measurement increases. This important property may be understood 
if one replaces the continuous fuzzy measurement by a sequence of instanta- 
neous fuzzy measurements with Gaussian weight functionals. It is plausible 
that such repetition of fuzzy measurements improves the resulting sharpness. 

The measurement of energy is now worked out as follows: Let the mea- 
sured system be initially in a pure state j^o)- Then for a particular given 

'In the preceding papers on RPI approach in continuous measurements an analogous 
parameter was usually called measurement error. However, use of this term may be mis- 
leading, and now we prefer a more definite term "resolution" . 
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readout [E] its further evolution is obtained by solving the corresponding 
Schrodinger equation (comp. (]TT[) with A = H and [a] = [E]) 

= [~H -k(h - E(t))^ (14) 

The norm of the final state vector \tpj<) gives then, according to the general 
formula (H), the probability density for the measurement readout [E\. 

P[E] = (VtI^t). (15) 

This together with the resulting state \ipT) is what can be known about the 
measurement. 

Expanding the state \ip t ) in the basis 

M*)> =e-* E » t / R |n> (16) 

of time-dependent eigenstates of H , we have the following system of equa- 
tions for the coefficients of the expansion \tp t ) = J2C n (t)\(p n (t)): 

C n = -K{E n - E(t)) 2 C n - l -J2(Vn\V\cp n >)C n ,. (17) 

n' 

The probability density of the measurement readouts is given by 

P[E}=J2\C n (T)\ 2 . (18) 

n 



3 "Free" multilevel system: model of deco- 
herence 



We shall consider first a "free" multilevel system with the Hamiltonian H Q . 
"Free" means vanishing V so that there is no driving field and no transitions 
between levels occur. We shall see in this simple case that a continuous fuzzy 
measurement of energy may serve as a model for decoherence leading to the 
resolution of the energy levels. 

Because of V — 0, Eq. (|17D have a simple solution 

((E n - Ef) T 



C n (T)= C n (0)exp 



— k / dt (E n - E(t)Y 



C n (0)exp 



AE 2 T 



(19) 



where ( }y is a time-average defined in (^) and AE? is defined by Eq. (pf) 
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3.1 Probability distribution 



A probabilistic analysis may be based on Eqs. 
probability density 

P[E]=J2\C n (T)\ 2 = E\ C nm 2 ew 



and (BD. In our case the 



-2 



({E n -Ef) T 



AE 2 T 



(20) 



is exponentially small when the temporal mean squared deflection ((E n — 
E) 2 )x of the function [E] from each level E n is more than AE 2 ^. On the 
contrary, the probability density is close to the maximum if ((E n — E) 2 )t is 
much smaller than AE 2 ^ for the level E n corresponding to a maximum value 
of|C n (0)| 2 . 

Let us now turn to the probability itself. The generalized unitarity has 
in the case of the energy measurement the form / d[E] W e JJ[e] = 1- Taking 
matrix elements of this equality in the energy representation, we have an 
equivalent condition 



d[E] exp 



,((E-E n ) 2 ) T 



AE 2 T 



1 (for all n). 



(21) 



The path integral here is of Gaussian type and can be evaluated. It may be 
shown that the (generalized) unitarity condition is fulfilled if the functional 



measure d[E\ is chosen to be d[E] = Y\ t y2Kdt/iTdE(t). 

Because of an exponential in the integrand, only paths close enough to 
the level E n contribute in ([H|). Namely, the integral which is analogous but 
restricted to the subset of readouts 



A n = {[E}\((E-E n ) 2 ) T < AE 2 T }, 



is close to unity: 



An 



d[E] exp 



,((E-E n ) 2 ) T 



AE 



r 



(22) 



(23) 



The probability that the measurement output [E] belongs to an arbitrary set 
A is according to the general formula (f|) and with the help of ([20]) 

Prob([£] e A) = [ d[E] Tr (U [E] p Uf E] ) = Y.Pn{A) |C n (0)| 2 (24) 
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with 

p n (A) = d[E] exp 

J A 

Comparing this with Eq. (|23|), we see that p n — 1 if A n C A. If, on the 
contrary, the set A does not intersect with A n (or their intersection is small 
compared with A n ) then p n = 0. 



-2 



((£-£n) 2 >T 



AE T 



(25) 



3.2 Different Regimes of Measurement 

As has already been said, the measurement resolution AEt decreases for 
increasing duration T of the measurement. To define different regimes of 
the measurement, this resolution should be compared with the characteristic 
difference AE ~ E n+ i — E n between energy levels.^ It is more convenient 
however to go over to the corresponding time parameters. Let T\ r be the 
duration at which the measurement resolution AEt achieves the value AE, 

T lr = — (26) 
kAE 2 K ' 

Then AEt is larger or smaller than AE depending on whether T is smaller 
or larger in comparison with T\ r . 

From Eq. (|19|) we have for the coefficients C n at time t = T\ r 



C n {T b ) = C n (0)exp 



{{E n -Ef) TlT 
AE 2 



(27) 



Therefore, the coefficient C n (t) decreases essentially during the time T\ r if the 
deflection of E(t) from E n is of the order or greater than AE. This is why 
the parameter T\ T may be called level resolution time. Additional arguments 
for this will be given below. 

Let us consider now different durations T and the corresponding regimes 
of measurement. In the case T ^> T\ r we have AEt "C AE. Now the 
probability density (|20| ) is not negligible only if for most of the measurement 
duration 

\E-E n \ < AEt<^AE (28) 

8 If the difference between levels substantially depends on the number n of the level, 
one can draw corresponding conclusions concerning the measurement in a certain energy 
band referring to the value AE which is typical for this band. 
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for some n. The corresponding set of readouts A n determined by Eq. ( p2|) is 
in this case a very narrow band close to the constant curve E(t) = E n . The 
width of this band is much less than AE, so that no other level E n > with 
n' ^ n lies inside A n . Taking the set of paths A in (|24j) in such a way that 
it includes A n but does not intersect with any other A n >, n' ^ n, we have 
p n = 1 and p n i = for n' ^ n. Therefore it follows from (p4|) that 



Prob(^l) = |C n (0)| 2 . (29) 

Thus, in the case T 3> T] r the measurement readout [E] will be an almost 
constant curve close to one of the levels E n , and the probability of it being 
close to a specified E n is |C n (0)| 2 . We see that the continuous fuzzy measure- 
ment, if it lasts longer than T\ r , distinguishes (resolves) between the levels. 
This justifies once more the name level resolution time for the parameter T\ T . 
To find the state of the system after the measurement in the considered case 
(T ^> Ti r ), we have to use Eq. ( [L9|) with the function E{t) satisfying the 



restriction (Pq). We obtain 

C n (T) = C n (0), C n ,(T) = 0, n'^n. (30) 

Thus, after a continuous fuzzy measurement of duration longer than T\ r 
the situation is similar to the one in the conventional description of the 
energy measurement as it is given in the von Neumann postulate: 1) The 
only possible measurement readouts [E] are those close to one of the energy 
levels E n , 2) the probability of [E] being close to a particular E n is |C n (0)| 2 
and 3) in the case of [E] being close to E n the system turns out to be in the 
energy eigenstate \(p n ) after the measurement. 

The physical process leading to the determination of an energy eigen- 
value (or of the eigenvalue of any other observable with discrete spectrum) 
has been investigated in the framework of different models and was called 
the process of decoherence [[T2|, |2T|. We see now that the continuous fuzzy 



measurement of energy fixes the energy level during T\ r . Therefore it presents 
a phenomenological description of the decoherence process leading to the de- 
termination of the specific energy eigenvalue. The characteristic mark of 
decoherence is the dying out the off-diagonal matrix elements of the density 
matrix. It is easily seen that it is this what really happens at the time T[ r 
because due to (p7| ) all products CiC* with % ^ j become exponentially small 
for T ^> T\ r . The time period during which non-diagonal elements die out, is 
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usually called the decoherence time. We call it level resolution time because 
we consider other regimes of measurement too, where the parameter T\ r plays 
another role. 

In this context the measurement performed during the time T <C T\ T might 
be considered as not yet finished. However this point of view is applicable 
only if the aim of the measurement is to distinguish between different levels. 
Continuous fuzzy measurement may however have as a goal to estimate the 
energy of the system with a resolution worse than AE. Then a period of 
measurement shorter than T\ T makes sense too. Let us turn to this regime. 

We consider the case T <C T\ T . This is equivalent to AEt 3> AE, so that 
the set of paths A n determined by Eq. (|22|) includes functions [E] in a very 
wide energy band around the constant curve E = E n . The width of this band 
is AEt, so that it includes many other energy levels. Vice versa, if we define 
the set of readouts A to contain all functions [E] in the energy band of the 
width AE T between -E m i n and E max , then many energy levels E n , n\ < n < n 2 
lie inside this band. For each of these levels intersection between A n and A 
makes a significant part of A n and therefore p n (A) ~ 1 for ni < n < n 2 . 
On the other hand, for n outside the interval [711,712] intersection between 
A n and A is small so that p n (A) is small. Making use of Eq. (^), we have 
approximately 

Prob(„4) = £ |C n (0)| 2 . (31) 

n=n\ 

This regime of measurement does not distinguish between the energy levels 
so that the character of the measurement in this regime is essentially the 
same as in the case of a continuous spectrum. 

4 The driven system: correlation between mea- 
surement readout and Rabi flopping 

We consider now a system which is exposed to an external driving field. For 
simplicity let it be a 2-level system under resonance influence of a periodical 
force with frequency uj = AE/h where AE = E 2 — E\. The potential 
energy has then the following non-zero matrix elements for the states (|16|): 
(<£>i|V|y> 2 ) = (v 9 2|V|v9 1 )* = V . Eq. ( fPTD gives for the related expansion 
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coefficients 



C x = -ivC 2 - k(E x - E(t)f d, 

C 2 = -ivC x - k{E 2 - E{t)f C 2 (32) 

with v = Vo/h. 

If no measurement takes place, i.e. for k = 0, such a system under- 
goes periodic transitions between the eigenstates (Rabi oscillations or Rabi 
flopping): 

C 1 (t)=R 1 (t), C 2 (t) = R 2 (t) (33) 

with 

R x (f) = d(0) cos vt-iC 2 (0) sin vt, 

R 2 (t) = C 2 (o) cosvt-id{0) siavt. (34) 

The parameter v = Vo/h is equal to a half of the Rabi frequency and the 
corresponding time scale is the Rabi period Tr = tt/v = hn/Vo- In the time 
Tr the system inverts completely from one level to the other and returns to 
the original level. 

Now the energy measurement is characterized by three different time 
scales: the level resolution time T lr , the Rabi period T R , and the duration 
of measurement T. With regard to physical interpretation one may look 
at them as follows: Small Tj r represents quick level resolution because of 
small fuzziness and therefore strong Zeno-type influence of the measurement. 
Small Tr corresponds to a strong influence of the driving field. For large T 
fuzziness effectively decreases. Let us consider different limiting cases for the 
relations between these times which will correspond to characteristic regimes 
of measurement. 

If T <C Tr, then the measurement duration is too short for the effect of 
the resonance influence to be observed. In this case we have in fact a "free" 
system as considered in Sect. |3|. The two possible relations, T\ r T and 
T T[ r , between the remaining parameters correspond to the cases when 
the measurement distinguishes between the levels and when it is too short 
to distinguish between themes discussed above in Sect. ||. 

We turn now to the opposite inequality, Tr < T. Now the influence 
of the driving field becomes stronger. But how big the resulting influence 
on the state of the system is, still depends on the amount of fuzziness of 
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the measurement. We will discuss in detail the regimes T\ r <C Tr (Zeno 
regime in which influence of the measurement is strong enough to damp out 
the Rabi oscillations), the opposite regime T R <C Tj r (Rabi regime in which 
Rabi oscillations maintain but there is no correlation between them and the 
measurement readout [E]) and the regime T R < 27rT lr but with T R , 2nT\ r and 
T of the same order. The latter is according to our goal the most interesting 
one, because the energy measurement readout turns out to be correlated to 
the oscillations of the state vector (modified Rabi oscillations). Some remarks 
concerning the regime T R C Ti r < T will be made. 



4.1 Zeno regime of measurement 

The first regime of measurement corresponds to the relation T\ T C Tr < T. 
One of its specific features results from the inequality T\ T <C T equivalent to 
AEt <C AE or kTAE 2 ^> 1. Because of this relation the damping terms 
in Eq. (|32]) act effectively when [E] differs from the corresponding level Ei 
or E 2 by a value of the order of AE or larger. This leads to a decay of 
the coefficients C\, C 2 . According to fllfif) the corresponding measurement 
readout [E] has low probability density. 

This decay of the coefficients may be prevented (and therefore the output 
may have large probability density) only if [E] is close to one of the levels 
Ei or E 2 . Recall our definition of being close. Intervals when E(t) deflects 
from the level more than by AE may occur because of the fuzziness of the 
measurement, but their complete duration should be much less than T lr and 
therefore much less than T. 

To verify that the probability density is actually large for and [E] close 
to E n , we consider a solution of Eq. ( |32l) for the case E{t) = E\.f\ Then 
Eq. (|3D reads: 

C x = -i^C 2 , C 2 = -i^C x -^C 2 . (35) 

J- R J R J lr 

These equations are easily solved to give 

d = a e Xlt + b e X2t , C 2 = 1— f Aia e Xlt + X 2 b e Mt ) (36) 



9 We simplify the consideration in this way, but the results are correct for the function 
[E] close, not identical to level E\ . 
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where 



TV 2 



Al ' 2 " ± \] 47j -J* (3?) 

and the coefficients a, b are determined by the initial conditions. For T\ T < 
Tr/27t we have exponentially decaying solutions. They can be completely 
described in a generic case, but we shall consider only the limiting regime 
Ti r <C Tr when the level resolution is much more rapid than the Rabi oscil- 
lation. We have in this case Ai = — 1/T\ T , A2 = —u 2 Xi, where v = itT^/Tr is 
a small parameter, and the coefficients take the form 



2 T 



d(T) = (Cx(0) -ivC 2 (0))e C 2 (T) = -ivC x {T). (38) 

Estimating the probability density by the formula flT8|), we find 



P(E = Ek) = (1Ci(0)| 2 + ^ 2 |C 2 (0)| 2 ) e~ 2v2 -r. (39) 

The exponential factor here is close to unity for a very wide range of T, up 
to the values of the order of T^/T\ T . Omitting this factor and neglecting the 
term containing z/ 2 , we have 

P[E = E 1 ] = \C 1 (0)\ 2 - (40) 
In the same approximation we have 

C l (T)=C l (0), C 2 {T) = 0. (41) 

We calculated the probability density for the constant curve (E = Ei) 
being the measurement readout. The same conclusion will be valid for the 
curves close but not identical to this constant curve. The set of such curves 
form a narrow band of the width AE?- Because of AE T AE this band 
does not include the second level E 2 . Just as in Sect. |372| , the probability 
of the measurement output to belong to this band is equal to | Ci (0) | 2 . It is 
clear that an analogous consideration may be applied for the measurement 
readout close to the energy level E 2 . 

Thus, for TJj C T R < T we have found a very simple picture: i) Only 
those measurement outputs [E] have high probability which are close to one 
of the constant curves E = E\ and E = E 2 \ ii) the probability of the output 
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to be close to E\ or E 2 is given by the initial values of the decomposition 
coefficients |Ci(0)| 2 or | C2 (0) | 2 correspondingly; iii) in case of the output 
being close to E\ or E 2 the final state is correspondingly the eigenstate \(pi) 
or \ip 2 )] iv) Rabi oscillations are completely damped away It is evident that 
this picture reflects the dominant influence of the quantum Zeno effect. It is 
therefore justified to call this regime the Zeno regime. 

4.2 Rabi regime of measurement 

We turn now to the case T R < T lr when Rabi oscillations are not damped. 
The three variants of this regime will be discussed. 

1. We consider first the case T R < T < T lr . What can be said about the 
measurement output [E]7 Returning to Eq. (|32|) , we have to apply now the 
inequality kTAE 2 <C 1 equivalent to T <C T\ r . It shows that the damping 
terms are essential only for very large deviation of [E] from both E\ and E 2 . 

If [E] lies in the band of the width AEt = AE^T\ r /T around E\, E 2 , no 
damping occurs. The solution of Eqs. (^) have then the form of the Rabi 
oscillations (|33| ) and the probability density for each of these measurement 
outputs [E] is equal (up to the order of magnitude) to unity: 

Prob[£] = |i?i(T)| 2 + \R 2 (T)\ 2 = 1. (42) 

All the curves [E] in the mentioned wide band may (approximately) be con- 
sidered as equally probable. Probability density essentially decreases only for 
the curves which deviate from both levels Ei, E 2 more than by the amount 
AEt (and therefore much more than by AE). The measurement in this 
regime does not respect the discrete character of the energy spectrum. The 
duration of the measurement is too short to resolve the levels. The energy is 
measured with the error AEt much greater than AE. 

Turning to the corresponding state of the system we note that with the 
damping terms being negligible in the present regime, Eqs. ( |52] ) have a very 
simple solution ([33|) describing Rabi oscillations. These oscillations are not 
prevented by the continuous measurement performed in this regime. They 
are sinusoidal with Rabi period Tr. We shall call therefore this regime of 
measurement the strong Rabi regime. The oscillations of the atom state are 
not reflected in this case in the readout [E]. 

2. The important variant of the Rabi regime is the one which shows 
a correlation between oscillations of the state vector and the measurement 
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readout. We shall demonstrate that this takes place for Tr < 27rTi r but with 
Tr, 27rT r and T being of the same order. . The analysis of Eqs. fl32|) is more 
complicated and less reliable in this case, but the process may be investigated 
with the help of the numerical simulation. 
Let us rewrite Eq. fl32| ) in the form 



e l {t) = K{E l -E{t)f = ^{^-^\ . (.11) 



C l = -ivC 2 -e 1 C 1 , C 2 = -ivd - e 2 C 2 (43) 

with 

f, (t) = k(E: - FAt^ 2 = 

T lr V AT 

Introduce instead of C n the new functions: 

S n (t) = e £ ^C n (t) (45) 

with 

£{t) = l -j\e 1 {t) + e 2 {t))dt. (46) 
Then these functions satisfy the equations 

S x + ivS 2 + -aSi = 0, S 2 + ivS 1 - \aS 2 = (47) 

whereby a = t\ — t 2 . Eqs. (|47|) imply the decoupled second-order equations: 

51 + {v 2 + \a - \a 2 ) S 1 = 0, 

5 2 + (v 2 - \a - la 2 ) S 2 = 0. (48) 

We shall assume (this will be justified a posteriori for the most probable 
measurement readout [T prob ]) that i) the deviation \E(t) — E\ with E = 
\{Ei + E 2 ) is not much larger than AT during most of the time, i.e. goes 
beyond these limits only for a short time or as a small deviation and ii) the 
spectrum of the function E(t) does not contain frequencies much higher than 
v (a half of the Rabi frequency). Then 

H < H < J^k_ (49) 
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Accordingly, to zeroth order in the ratio Th/2ttT\ t we neglect the two last 



terms in the brackets of Eq. (48). The functions S n are then approximated 
by harmonic functions of frequency v, i.e. by the Rabi oscillations R n of (51). 
We have therefore zeroth approximation 

C x {t) = e- e ®R x (t\ C 2 (t) = e- £ ^R 2 (t). (50) 

This approximation reflects the general oscillatory character of the evolution 
of the state \ipt), but it is too rough to use it for calculating the probability 
and for estimating the most probable measurement readout [E] . Instead, we 
shall use another procedure for this. 

Consider again the functions C n (t) at an arbitrary time moment < t < 
T and form the function 

P(t) = \C 1 (t)\ 2 + \C 2 (t)\ 2 . (51) 

The quantity P(t) depends on the values {E(t')\0 < t' < t} of the function 
[E] in the moments preceding t. This quantity can be interpreted as the prob- 
ability density of the corresponding measurement readout if the measurement 
is finished at time t. With increasing time t the value P(t) decreases. 

Let us analyze how the function P(t) depends on time and choose the 
function [E] in such a way that P(t) decreases as slowly as possible. Then 
the resulting function [E] will correspond to the most probable measurement 
readout. 

Differentiating Eq. ([y]) and making use of the Eq. (|32l ) we have 



P(t) = -2k - E(t)f |d| 2 + (E 2 - E(t)f \C 2 \ 2 \ . (52) 

We minimize the absolute value of the r.h.s. and obtain for the most probable 
measurement output 



E 1 \C 1 {t)\ 2 - 


-E 2 \C 2 (t)\ 2 




-|c 2 (t)| 2 



^ob(t) = / , 2 ; v^r^r ■ (53) 



Now we can replace C n by the approximate solution (|50|) . This gives for 
the most probable measurement output 

E pmh (t) = E 1 \R 1 (t)\ 2 + E 2 \R 2 (t)\ 2 . (54) 
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This curve is correlated with the Rabi oscillations of the state in the sense 
that the curve E pro b(t) is closer to the level E n when the probability for the 
system to be on this level is larger. The function E proh (t) oscillates around 
the middle line E between the levels. The amplitude of the oscillations 
is equal to |A£y '((C^O)! 2 - |C 2 (0)| 2 ) 2 + 4(Im(CiC 2 *)) 2 . If the initial value 
of the decomposition coefficient Ci(0) is equal to 6*2(0), the most probable 
measurement output -E pro b coincides with the constant curve E = E. Only 
in this case the Rabi oscillations disappear. 

The measurement readout (|54]) is shown to be most probable. However, 
it is not clear from the above argument how rapidly the probability density 
decreases if the readout [E] deviates from its most probable shape. If it 
decreases rapidly enough, then we may speak about correlation between the 
state oscillations and measurement readout. This is investigated with the 
help of the numerical calculation. We omit details. The following conclusions 
are obtained: 

i) Oscillations of the state vector between the energy eigenstates (Rabi regime) 

takes place for Tr < 2ttTi t . The boundary to the Zeno regime is marked 
by T R « 2nT h . 

ii) The state oscillations are modified Rabi oscillations. Their frequency 

becomes smaller than the Rabi frequency when Tr approaches 2-kT^. 
Together with this the pure sinusoidal shape of the oscillations is lost. 
This indicates that there is already some Zeno-type influence of the 
measurement. See the figure for more details. 

iii) Correlation between the modified Rabi oscillations and the measurement 

readout [E] is substantial for Tr close to 2ttTi t . In the other case the 
band of different readouts [E] having large probability becomes too 
wide to speak about a correlation. 

3. The complex case Tr C Ti r < T is not treated in detail here. We 
restrict ourselves to some comments. In this case small errors related to the 
approximate character of (|50"D grow exponentially, so that this solution is 
inappropriate for the analysis, and the conclusions made above are not valid. 
The evolution is qualitatively different in this case. Under these conditions 
the system "forgets" because of the measurement its initial conditions and 
the oscillating curve ([53]) cannot arise as a typical measurement output. 
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The mechanism of "forgetting" is loss of coherence leading to a mixture 
instead of a pure state. If the duration T of the measurement is not too 
long, the result of the measurement and the evolution of the system may 
be adequately presented by a single readout [E] and the resulting pure final 
state \ip [ T ] ). A single readout is a set of null measure and strictly speaking 
is characterized by zero probability. One has to consider a band of readouts 
instead and integrate the probability density P[E] over this band to obtain 
non-zero probability. However close measurement readouts [E] lead to close 
final states of the system if the measurement is not too long. Therefore, the 
final state corresponding to the band of readouts will be the same pure state 

\E] 

\ipT ) which corresponds to each of the readouts in the band. 

This argument is not valid however for a long measurement. In this case 
even very close readouts [E] lead to different final states IV't'')- Therefore, 
even a very thin band of readouts gives a mixed final state. With increasing 
T this mixed state becomes universal, not depending on the initial state of 
the system. 

Indirectly the latter statement is confirmed by the non-selective descrip- 
tion of the measurement with the help of the density matrix (BJ). It has been 



shown in |22| that this density matrix satisfies a master equation having the 
form of the Liouville equation but with an additional term proportional to 
the double commutator [AL4p]]. In our case this equation reads as follows 

P=-±[H q + V,p]-%[Ho[Hqp]] (55) 

and can easily be solved for a resonantly driven two-level system. This gives 
for 87rTi r > Tr (Rabi regime) 

pn - P22 = e~ ( ' /4Tlr) (acosfit + b sin fit). (56) 



Here matrix elements are taken between the states (16) , = y (2v ) 2 — (2Ti r )~ 2 
and the parameters a, b are determined by initial conditions. It is evident 
that the difference p\\ — P22 due to non-symmetrical initial conditions disap- 
pears during time of the order of T\ r . 
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5 Possible realizations of continuous fuzzy mea- 
surements of energy 

We want to sketch briefly possible realizations of continuous fuzzy measure- 
ments of energy. Our suggestions need more detailed elaboration which 
cannot be given here. An output of an instantaneous measurement is a 
number. An output of a sequential measurement is a series of numbers 
ax, a,2, ■ ■ ■ , a n , .... If the time scale is larger than the interval between in- 
stantaneous measurements this series of numbers may be identified with a 
continuous curve a(t) presenting the readout of a continuous measurement 
(monitoring) . Fuzziness of the continuous measurement depends on the fuzzi- 
ness of the instantaneous measurements and on the intervals between them. 
We give three examples: 

1. The most obvious approach is to make a continuous sharp measurement 
fuzzy in making the underlying instantaneous measurements fuzzy, in other 
words, to modify measurements leading to the quantum Zeno effect in this 



way. An arrangement for the Zeno effect was suggested in [23] and used in [0]. 
A system (ion in a trap) consisting of levels 1 and 2 is driven between these 
levels by a resonant perturbation. There is a subsidiary level 3 that can decay 
only to level 1. The state measurement is carried out by action of an optical 
pulse in resonance between levels 1 and 3. This pulse can induce the 1 — > 3 
transition with the subsequent spontaneous return to level 1 accompanied 
by the emission of a photon with the same frequency. It causes a projection 
of the ion onto level 1 accompanied by scattering the photon. This happens 
with the probability |ci| 2 if the initial state of the ion is ci|l) + C2|2). With 
the probability | C2 1 2 the pulse projects the ion onto level 2 with no scattering. 

The measurement of this type is sharp (of von Neumann type) and an- 
swers the question whether the ion is on level 1. Frequent repetition of this 
measurement was shown to lead to the quantum Zeno effect freezing the 
atom on level 1. Each measurement in a series completely resolves in this 
case between the energy levels, i.e. is not fuzzy. It is necessary to modify 
the experimental setup to realize a fuzzy measurement of energy.0 

To do this, one may choose an ion with very fast decay of the subsidiary 

10 The repeated measurement of the same type but with varying parameters of the system 
has been considered in [ [To| . It is however difficult to compare the results obtained in this 
paper with our results because fuzziness was not the aim of investigation in ]lQ] . 
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level 3 to both main levels 1 and 2 and a continuous spectrum of the radia- 
tion inducing the transitions 1 — >• 3 and 2 — > 3. The width of level 3 and the 
width of the spectrum of the inducing radiation have to be more than the 
difference AE between levels 1 and 2. The induced transition from levels 1 
and 2 to the subsidiary level 3 with the subsequent spontaneous transition 
back to the same level will be accompanied by elastic scattering of a photon. 
Measuring the frequency of the scattered photon is equivalent in these con- 
ditions to an instantaneous fuzzy measurement of energy (performed with 
the resolution worse than AE). A series of scatterings is then a continuous 
fuzzy measurement of energy. Resolution of the continuous measurement 
AEt depends on the length of the series or duration T of the continuous 
measurement. Special precaution should be taken in this scheme to provide 
that only transitions 1 — > 3 — > 1 and 2 — > 3 — > 2 might occur. For this, level 
3 may be double, consisting of two very close levels V and 2' with the only 
transitions 1 V and 2 <->• 2' permitted. 

2. An electron in a magnetic field is a two-level system and repeated 
measurements of its energy can be performed with the help of several Stern- 
Gerlach (SG) devices. The position of the electron after scattering in the SG 
device depends on its energy. Therefore measuring the final position of the 
electron is in fact measuring its energy. If the initial position is uncertain 
to some degree, the resulting measurement of energy is fuzzy. In fixing the 
initial uncertainty, we can determine the fuzziness of the instantaneous mea- 
surement of energy. If now the electron is scattered subsequently in a series 
of SG devices, we have in fact a series of fuzzy instantaneous measurements 
of energy approximating a continuous one. 

As a driving field we can include an external electromagnetic field in res- 
onance with the difference between the electron levels. This induces periodic 
transitions between the levels. Then all the effects discussed in Sect. f| can 
be observed with such a system. An analysis is given in ||, |!|. The results 
are in agreement with our general results. 

3. One more way of introducing fuzziness can be found in quantum 
optics. The multilevel "free" system is in this case realized as a mode of the 
electromagnetic field. The measurement of energy of this system (usually 
referred to as the measurement of the photon number) is realized by the 
interaction of this mode with another (probe) electromagnetic wave through 



a nonlinear optical Kerr medium in which both waves propagate [ft4 |. The 



24 



interaction Hamiltonian has the form 



Hj = h x (a+a) (6+6) (57) 

where a, 6 are the annihilation operators of the two modes. The observable 
to be measured indirectly is H = a^a, but the directly measured observable 
is a quadrature phase component 

6 2 = i(6-6t). (58) 

After the probe wave has travelled the distance L in the Kerr medium, 
the value of this component is, due to the interaction Hj, 

6 2 (L) ~ ^(6 1 (0))// + A6 2 (0) (59) 

where A6 2 (0) is an initial uncertainty of the amplitude 6 2 and (6i(0)) an 
initial mean value of the second quadrature component b\ = \{b + 6+). This 
means that the output for the measurement of H is read off as the value 
of the component 6 2 . Fuzziness of the measurement depends therefore on 
A6 2 (0) and can be controlled. The series of the measurements of this type 
may be considered as a continuous fuzzy measurement of energy H . 



6 Conclusions 

Fuzziness of a measurement is caused by the low resolution of the detec- 
tor. Such a measurement is unable to transfer the state of the system to an 
eigenstate of the measured observable, and the measurement readout must 
not agree with an eigenvalue. Such a continuous fuzzy measurement as mon- 
itoring an observable may be approximated by as a series of unsharp or 
imperfect "instantaneous" measurements. It is qualitatively different from a 
series of sharp (von Neumann type) measurements. 

We introduced fuzziness with the help of restricted path integrals, on 
the basis of the resulting non-Hermitian Hamiltonians. We discussed the 
continuous fuzzy measurement (monitoring) of energy for a multilevel and in 
more detail for a two-level system leading to a readout E{t). In addition to 
the duration T of the measurement there are two other characteristic time 
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scales. The level resolution time T\ r represents the strength of an inherent 
Zeno-type influence ("level freezing"). If an external resonant driving field 
is acting upon the system, the strength of its influence is represented by 
the Rabi period T R . According to the relative orders of magnitudes of these 
times, different regimes with characteristic properties of the measurement 
readout and the related evolution of the state vector can be distinguished. 
We obtained the following results: 

i) The study of a "free" multilevel system which is only under the influence 

of the measurement (no driving field present) shows that continuous 
fuzzy measurement of energy is a model for decoherence leading to the 
resolution of separate levels. 

ii) The discussion of a resonantly driven two-level system reveals two dif- 

ferent regimes of measurement: In the Zeno regime the oscillations of 
the state vector are prevented and the measurement readouts are es- 
sentially discrete. In the Rabi regime the state oscillations maintain 
and the readouts may be arbitrary curves. 

iii) An important new effect is predicted: For time scales corresponding to 

the Rabi regime but being close to the border with the Zeno regime 
the frequency of the state oscillations is decreased and their shape de- 
viates from pure Rabi oscillations. Correlation are found between the 
oscillations of the state vector and the measurement readouts. Thus, 
in following E(t) one can directly read off the oscillations of the state 
vector. 

Possible ways of realization of a fuzzy continuous measurement of energy 
are shortly commented. 
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Figure capture 

Correlation between modified Rabi oscillations of the state of a resonantly 
driven two-level system and the readouts of a continuous fuzzy measurement 
of energy. We assume T R /27rT lr = 0.8. a) The in the upper diagram shows the 
periodic oscillations of the state vector between the two energy eigenstates. 
The period is close to but larger than the Rabi period Tr. The shape is of 
a slightly modified Rabi oscillation, b) The diagram in the middle shows as 
the dashed line the most probable energy readout [-Eprob] which is strongly 
correlated to the state oscillation. Two bands around it (with the width 
W = 2 and W = 3 are drawn. c)The probability P that an energy readout 
[E] lies in the band with the width W around the most probable readout is 
given in the lower diagram. 
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